Universal Measure of Entanglement 
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A general framework is developed for separating classical and quantum correlations in a multi- 
partite system. Entanglement is defined as the difference in the correlation information encoded by 
the state of a system and a suitably defined separable state with the same marginals. A general- 
ization of the Schmidt decomposition is developed to implement the separation of correlations for 
any pure, multipartite state. The measure based on this decomposition is a generalization of the 
entanglement of formation to multipartite systems, provides an upper bound for the relative entropy 
of entanglement, and is directly computable on pure states. The example of pure three-qubit states 
is analyzed in detail, and a classification based on minimal, four-term decompositions is developed. 
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Entanglement is the simultaneous occurrence of super- 
position and correlation in composite systems, and gives 
rise to some of the most counterintuitive phenomena of 
quantum mechanics. It was Einstein's early misgivings 
about the nonlocality implicit in entangled states jlj that 
led to the celebrated EPR paper j^]. This paper in 
turn inspired Schrodinger's cat paradox and his prophetic 
identification of entanglement as the most characteristic 
feature of quantum mechanics A striking confirma- 
tion of this prediction is provided by Bell's renowned in- 
equality, which is violated by a pure state if and only if 
it is entangled 0. Long at the center of most founda- 
tional discussions, and widely regarded as the essence of 
quantum strangeness, entanglement has in recent years 
emerged as the key resource in quantum information pro- 
cessing 1^. This remarkable transformation of paradox 
into paradigm was brought about by the explosion of 
interest in quantum computing and communication in 
recent years 0. 

The quantification of entanglement has thus emerged 
as a central problem in quantum information theory. 
Among the proposed measures, the pioneering contribu- 
tions of Bennett et al. 0, Q defined entanglements of 
formation and distillation, Ep and Eo respectively, on 
considerations of convertibility vis-a-vis maximally en- 
tangled pairs. Indeed Ep and E]j (< Ep) are in effect 
benchmark buy and sell exchange rates for maximally 
entangled states, with any other "good" measure of en- 
tanglement expected to fall in between 0. On the other 
hand, the relative entropy of entanglement To'l, denoted 
by Eh, is based on distinguishability from the set of sepa- 
rable states, with relative entropy serving as a measure of 
contrast |ll| . These and most other proposed measures 
are based on operational considerations of quantum com- 
munications or similarly motivated mathematical axioms 
[l^ . While bipartite entanglement is for the most part 
well understood, a general formulation of the multipar- 
tite case is lacking and remains an outstanding unsolved 
problem. 

In this Letter we develop a definition of entanglement 



based on general, rather than operational or axiomatic, 
considerations. We also develop a generalization of the 
Schmidt decomposition [l^ to implement this separation 
for pure, multipartite states, and extend this to all multi- 
partite states. To see the main idea in qualitative terms, 
let us consider a state p of a composite system together 
with its reduced (i.e., unipartite) density matrices p^, 
, . . p^ , which represent the states of the components. 
The latter are the quantum marginals of p, and collec- 
tively carry less information than does p, the difference 
being the correlation information encoded by p. Further- 
more, because of the superpositions it may contain, p in 
general encodes more correlation information than any 
classically correlated state that possesses the same set 
of marginals. The difference, minimized subject to suit- 
able constraints, is then a measure of nonclassical corre- 
lations, or entanglement, encoded by p. Note that this 
definition is based on separating the correlations encoded 
by a quantum state into classical and nonclassical parts 
by reference to a common set of marginals |l4l| . Note also 
that this notion of separation is already implicit in the 
standard entanglement measure for pure bipartite states. 

We now proceed to a quantitative formulation of the 
above ideas. Our objective is to define a classically cor- 
related state a associated with p. We require a to (i) 
be separable, i.e., a mixture of pure, product states, and 
(ii) possess the same set of marginals as p. We further 
require a to encode the greatest correlation information 
while having the least contrast ^3 against p. Constraint 
(i) ensures that the correlations encoded by a are classical 
in nature , while constraint (ii) guarantees that these 
correlations can be meaningfully compared to those of p. 
Furthermore, the maximum correlation condition stipu- 
lates that a capture as much of the correlation content 
of p as possible, so that only nonclassical correlations are 
counted as entanglement, while the minimum contrast re- 
quirement requires that a differ as little as possible from 
p. Following Ref. |To| . we use relative entropy as our 
measure of contrast. 

For pure bipartite states, minimizing relative entropy 
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subject to constraint (i) is sufficient to satisfy (ii) and 
the maximum correlation condition as well [lOj. How- 
ever, the bipartite case is special, as evidenced by, e.g., 
the fact that the Schmidt decomposition does not 
in general extend beyond N = 2 [ig. The distinguish- 
ing characteristic here is the index Nms, which we de- 
fine as the number of distinct marginal spectra possessed 
by a pure state Now Nms = 1 for N = 2, and 

1 < Nms < N ior N > 3. Furthermore, one can easily 
show that a pure state is Schmidt-decomposable if and 
only if Nms = 1 ■ Thus A'' = 2 is special, and there 
is no a priori reason to expect the maximum correlation 
and minimum contrast requirements to be equivalent in 
the multipartite case. Since this equivalency is a pre- 
requisite for the above optimization problem, we add the 
condition tr[(cr — p) log(cr)] = 0, labelled (iii), which turns 
out to have an interesting interpretation as seen below. 

Let M{p) be the set of states that satisfy constraints 
(i)-(iii). Then the correlation information contained in p 
and a e M{p) are respectively given by C = S{p) — S{p) 
and C — S{a) — S{a), where S{x) stands for the von 

Neumann entropy of x, and p=a p^ ® p^ ® . . . ® p^ 
represents the joint, uncorrelated state of the marginals 
of p or cr. The nonclassical correlation information is 
thus given by the difference CP - C = S(a) ~ S{p). 
Since the above extremal conditions require the simulta- 
neous maximization of C"^ and minimization of S{p\\a), 
we must ensure that their sum, S{d) — 8(0') + S{p\\cr), is a 
fixed quantity. This is precisely what is achieved by (iii), 
which reduces the sum to S{p) — S{p), thus guaranteeing 
the desired equivalency. Moreover, for cr G M{p), (iii) is 
equivalent to 



S{p\W) + S{a\\p) = S{p\\p), 



(1) 



which stipulates an additivity condition on contrast in 
comparing p to p via a. Note that inasmuch as contrast 
{— relative entropy) plays the role of distance here, we 
may interpret (iii) as requiring that ct be on a straight 
contrast line from p to p. 

The entanglement measure is now defined as the min- 
imum of nonclassical correlations: 



E{p) mm, 



aeM{p) [S{cr) - Sip)] 
min^e^(p) S{p\\ a), 



(2) 



and will simply be referred to as "entanglement." Note 
that because of (iii), A4{p) is not convex, so that there 
may be multiple local minima for (0). Note also that 
E = Ep = Eji on pure bipartite states, and because 
of (ii) and (iii), E > Eji on all states. This completes 
our discussion of the general framework and definition of 
entanglement. 

Next, we develop a generalization of the Schmidt de- 
composition [l9| as a means of developing an entangle- 
ment measure on pure, multipartite states. Specifically, 



our objective is a minimal, orthogonal decomposition of 
any pure state lip) such that its decohered form (t('0) 
obeys conditions (i)-(iu); we shall refer to these as com- 
pact decompositions. Compact decompositions are in ef- 
fect continued Schmidt decompositions: start by apply- 
ing the Schmidt decomposition to Itp) with respect to A 
and (BC . . . Z) as the two subsystems. 



1^) = E, 



l,BC...Z\ 



(3) 



then to l^j^*-^ "^) with respect to B and {C . . . Z) as the 
two subsystems, and continue the process to the last step 
where Y and Z are decomposed. The result is 



A \B 



YZ N 



Y 



(4) 



This expansion has a tree structure with each ket and 
its associated coefficient carrying a set of indices which 
label the successive branch numbers. Accordingly, the 
collection of kets with all but the rightmost index (which 
is set off by a semicolon) equal constitute an orthonormal 
set. For example. 
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= 5. 



(5) 



For each such condition, there is a corresponding one for 
its Schmidt mate, namely 
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XYZ 



XYZ 



= 6. 



(6) 



Similarly, the collection of A's with all but the rightmost 
index fixed constitutes a probability set. For example. 



A, 



w 



1. 



(7) 



Note that the last two subsystems in Q share a com- 
mon set of indices and a common coefficient as in a stan- 
dard Schmidt decomposition. Thus there are at most 
d-^ d^ . . . d-^ mm{d^ , d^ ) terms in the above expansion, 
where, e.g., d"^ is the dimension ofH^, the Hilbert space 
belonging to subsystem A. Except for the last pair, 
changing the order of subsystems in Eq. Q in general 
results in a different decomposition. Therefore, allowing 
for possible degeneracies, there are N\/2 compact decom- 
positions for an A^-partite system. For three qubits, e.g., 
there are no more than three decompositions each con- 
taining no more than 4 terms. 

The decohered form of -0 in the basis of Eq. Q is 



^ — ^7, 



A \B 



Af A 



.A 



YZ 



^t^b'»■■■'»^l^b■■■^.■,^y 



(8) 



where each state vector in Q has been replaced with 
the corresponding projection operator in ijSJ). The or- 
thonormality conditions stated above can now be used to 
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show that (a) the expansion given in Eq. © is orthog- 
onal, (b) a^-^'"^{ip) has the same marginals as tp, and 
(c) tr[CT^^-^(V') - |V')(V'|]log[cr'^^-^(V')] = 0. In other 
words, a^^---^{ip) satisfies conditions (i)-(iii) above. It 
is not difhcult to show that among all possible ways 
of decomposing a pure state by means of a sequence 
of Schmidt decompositions, compact decompositions are 
unique in guaranteeing conditions (a)-(c) stated above. 

Denoting the set of N\/2 possible cr's by A^^(V'), we 
define the entanglement of ip following Eq. 



where 16*+) =^ cos(6l/2)|+) + sin(6l/2)|-) and 16* 
sin(6'/2)|-|-) — cos(6'/2)|— ). In addition to the normal 
ization condition ^^Pi — 1, there are two further con- 
straints on the four real angles appearing in 114|l . 



dof 



dcf 



(9) 



(10) 



min^g^c(^) 5(cr), 

where the superscript "c" refers to compact decomposi- 
tions as the basis of E'^. If this minimum is realized on 
say, a"^^'"^ (tp), we have 

log(Af„Af^^,,...Ari. ^ 

The bounds on the entanglement of a system of iV par- 
ties each with a Hilbert space of dimension d are given 
hy Q< E" <{N -I) log(d). Note that Eqs. © and fTn| 
provide a readily computable, finite prescription for find- 
ing E'^i^ip). This is in contrast to Eqs. (O and ifTTjl which 
involve nontrivial minimizations over infinite sets. It can 
also be verified that E'^ (as well as E) is invariant under 
local unitary operations, and that its expected value is 
non-increasing under local operations and classical com- 
munications 8, 10]. 

Eq. (|10|l can be rearranged as 



E^{^)^Hi\^)+J2_ \tM\l) 



A \B 



A. A, 



H{X1...J, (11) 



where, e.g.. 



dof 



Af 



log(A: 



(12) 



is the information associated with the probability set 
{Af .}. This nested pattern is of course inherited from the 
tree structure of compact decompositions noted above. 

Next we consider the thrce-qubit case as an illustration 
of the foregoing construction. There are three compact 
decompositions for a tripartite pure state ip which may 
be labelled ip^^'^ , ip^'^^, and ip'^^^ , where, e.g., 



ABC\ 



(13) 

For three qubits represented as spin states, each of these 
can be brought to the following standard form by local 
unitary transformations: 



kABC 



Pi 
P2 e 
Pa 

pl e 



if3\ \A 



tan(6i{,/2) tan(6lc/2) 



jPiPs)^ + (P2P4)^e'(^-") 
(piP4)^e*'5 -I- {p2P?.)^e~ 



(15) 

leaving a total of five independent, real parameters. This 
is in agreement with previous results [l9j , although these 
involve a minimum of five terms in their decompositions 
as compared to four here. 

When decohered, each of the three compact decompo- 
sitions of i) yields a diagonal, separable density matrix 
of no more than four terms which satisfies conditions (i)- 
(iii). Among these, the density matrix with the least 
entropy defines the entanglement of i>. Therefore, E'^{ip) 
equals —J2iPi log(K)' O'' simply H{p*), where {p*} is the 
probability set corresponding to the minimal state. It is 
clear from this formula that < E'^{ij}) < log(4). A sim- 
ple example of a maximally entangled three-qubit state 
is obtained from (|14|l by setting a — 0, /3 — tt, db — 0, 
dc ~ TT, and p* = 1/4. The three decompositions of this 
state are identical due to its symmetry. This state may 
be represented as 



®i(| 
®i(|- 



BC 



BC 



(16) 



\-r ® \eb-r ® \e,-r , (w) 



The entanglement of this state is log(4) units which, by 
way of comparison, is twice as large as the entanglement 

of the GHZ state i(| + ++) + \ )) H^. Referring to 

Eq. Hll|) . we interpret this result as representing log(2) 
units for the bipartite entanglement of A versus (BC), 
and log(2) units for the average entanglement of the two 
Bell states that make up p^'~^ . 

Pure three-qubit states are naturally classified in terms 
of the ranks of their marginal spectra (see Ref . [l^l for an 
alternative scheme) . Class I corresponds to all marginals 
being of rank one, and consists of product states with zero 
entanglement. These correspond to = for z 7^ 1 in 
(|14|l . with Nms — 1. Class II corresponds to one marginal 
being of rank one, the other two full rank. These are 
product states of pure one-qubit and two-qubit states, 
with the entanglement carried by the latter. Such states 
correspond to ps = p4 — in ifHI) . with N^s — 2. 
Class III corresponds to all marginals being full rank. 
This class is conveniently subclassified as III-a,b,c, corre- 
sponding to Nms = 1,2,3, respectively. With N^s = 1, 
Class Ill-a admits a standard Schmidt decomposition [l^ 
corresponding to p2 — P3 = 9b = Oc = in (fT^ . The 
further restriction pi — p4 leads to the GHZ state men- 
tioned above. Class Ill-b has Nms = 2 and corresponds 
to 6b = Oc in H14|) . Class III-c with Nms = 3 is the typical 
case shown in (I14II . It is clear that only class III states 
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carry tripartite entanglement, while I and II are better 
described as products of one-qubit and two-qubit states. 

Having described compact decompositions and the as- 
sociated entanglement measure E'^ on pure, multipartite 
states, we now proceed to extend the latter to mixed 
states. The method we use was pioneered in Ref. 0, 
and has since been found to be a generally valid strategy 
for extending pure state measures to all states The 
key idea is to define the entanglement of an ensemble of 
pure states to be the ensemble average of the entangle- 
ment. The minimum of this quantity over all possible 
ensemble representations of a mixed state is then defined 
as its entanglement. In symbols, 

E^ip) minV (17) 

^ — 

where the minimum is taken over all possible ensemble 
representations X]a^'alV'a)(^a| of p. It is interesting to 
note that ensemble averages of pure state entanglement 
are implicit in the tree structure of compact decomposi- 
tions, as may be seen in Eq. and in the above ex- 
ample of a maximally entangled three-qubit state as well. 
Note also that since E'^ and Ep are equal on pure bi- 
partite states, they must coincide on all bipartite states. 
Therefore, E'^ as extended by Eq. (|17|l is a generalization 
of Ep to all multipartite states. Furthermore, because 
compact decompositions satisfy conditions (ii) and (iii), 
we have E'^ > En on all multipartite states. 

Thus far we have defined two measures of entangle- 
ment, E in Eq. |(2Jl, and in Eqs. © and lO. Consid- 
ering pure states for the moment, we defined these mea- 
sures by minimizing nonclassical correlations over and 
TW^, sets of states that satisfy conditions (i)-(iii). While 
Ai includes all such states and is infinite, Ai'^ includes 
states derived from compact decompositions only and is 
finite. Thus E < E'^ on pure states, with equality known 
to hold on bipartite states. Since we expect compact de- 
compositions to be minimal representations, it is likely 
that E'^ — E ox\ all pure states, although we have not 
succeeded in establishing this equality in general. In the 
event that the equality holds on pure states, there arises 
the question of whether the two measures coincide on 
mixed states as well, also an open question at present. If 
E turns out to be strictly smaller than E'^ on pure states, 
then E'^ provides a readily computable upper bound for 
E. 

Entanglement has been defined here as nonclassical 
correlation information, without any direct reference to 
operational considerations "2^. The results, however, are 
remarkably close to operationally inspired measures of 
entanglement, affirming the expectation that there must 
be a universal meaning to the quantity of entanglement 
embodied in a quantum system. 

This work was supported in part by a research grant 
from California State University, Sacramento. 
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